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Abstract 

Every three-connected planar graph with n vertices has a drawing on an 0{n^) x O(n^) 
grid in which all faces are strictly convex polygons. These drawings are obtained by 
perturbing (not strictly) convex drawings on 0(n) x 0(n) grids. Tighter bounds are 
obtained when the faces have fewer sides. In the proof, we derive an explicit lower bound 
on the number of primitive vectors in a triangle. 

1 Introduction 

A strictly convex drawing of a planar graph is a drawing with straight edges in which all faces, 
including the outer face, are strictly convex polygons, i. e., polygons whose interior angles are 
less than 180°. 

Theorem 1. (i) A three- connected planar graph with n vertices in which every face has 
at most k edges has a strictly convex drawing on an 0{nw) x 0{n^k/w) grid of area 
0{n^k), for any choice of a parameter w in the range 1 < w < k. 

(ii) In particular, every three- connected planar graph with n vertices has a strictly convex 
drawing on an 0{n^) x O(n^) grid, and on an 0{n) x 0{n^) grid. 

(iii) For k < 4, an 0{n) x 0{n) grid suffices. 
The drawings can be constructed in linear time. 

When referring to a W x H grid of width W and height H, the constant hidden in the 
O-notation is on the order of 100 for the width and on the order of 10000 for the height. This 
is far too much for applications where one wants to draw graphs on a computer screen, for 
example. For the case w = 1, the bound is tighter: the grid size is approximately lAn x 30n^k. 
For part (iii) of the theorem, the grid size is at most 14n x 14n, and if the outer face is a 
triangle, it is 2n x 2n. 

The main idea of the proof is to start with a (non-strictly) convex embedding, in which 
angles of 180° are allowed, and to perturb the vertices to obtain strict convexity. We will use 
an embedding with special properties that is provided by the so-called Schnyder embeddings, 
which are introduced in Sectional 

*Renyi Institute, Hungarian Academy of Sciences, PoB 127, Budapest 1364, Hungary, barany@renyi.hu, 
and Mathematics, University College London, Gower Street, London WCIE 6BT, UK. 

^Freie Universitat Berlin, Institut fiir Informatik, Takustrafie 9, 14195 Berlin, Germany, rote@inf.fu- 
berlin.de 



Historic context. The problem of drawing graphs with straight hnes has a long history. 
It is related to realizing three-connected planar graphs as three-dimensional polyhedra. By 
a suitable projection on a plane, one obtains from a polyhedron a straight-line drawing, a 
so-called Schlegel diagram. The faces in such a drawing are automatically strictly convex. By 
a projective transformation, it can be arranged that the projection along a coordinate axis 
is possible, and hence a suitable realization as a grid polytope gives rise to a grid drawing 
of the graph. However, the problem of realizing a graph as a polytope is more restricted: 
not every drawing with strictly convex faces is the projection of a polytope. In fact, there is 
an exponential gap between the known grid size for strictly convex planar drawings and for 
polytopes in space. 

The approaches for realizing a graph as a polytope or for drawing it in the plane come 
in several flavors. The classical methods of Steinitz (for polytopes) and Fary and Wagner 
(for graphs) work incrementally, making local modifications to the graph and adapting the 
geometric structure accordingly. Tutte ^3 El gave a "one-shot" approach for drawing graphs 
that sets up a system of equations. This method yields also a polytope via the Maxwell- 
Cremona correspondence, see jTTj . All these methods give embeddings that can be drawn on 
an integer grid but require an exponential grid size (or even larger, if one is not careful). 

The first methods for straight-line drawings of graphs on an 0{n) x 0{n) grid were pro- 
posed for triangulated graphs, independently by de Fraysseix, Pach and Pollack and by 
Schnyder JSI- The method of de Fraysseix, Pach and Pollack ffl is incremental: it inserts 
vertices in a special order, and modifies a partial grid drawing to accommodate new vertices. 
In contrast, Schnyder's method is another "one-shot" method: it constructs some combina- 
torial structure in the graph, from which the coordinates of the embedding can be readily 
determined afterwards. Both methods work in linear time. 0{n) x 0{n) is still the best 
known asymptotic bound on the size of planar grid drawings. 

If graphs are not triangulated, the first challenge is to get faces which are convex. (Without 
the convexity requirement one can just add edges until the graph becomes triangulated, draw 
the triangulated supergraph and remove the extra edges from the drawing.) Many algorithms 
are now known that construct convex (but not necessarily strictly convex) drawings with 
0{n) X 0(n) size, for example by Chrobak and Kant [3] (a la Fraysseix, Pach and Pollack); or 
Schnyder and Trotter J3j and Felsner 0, see also (a la Schnyder). Our algorithm builds 
on the output of Felsner's algorithm, which is described in the next section. Luckily, this 
embedding has some special features, which our algorithm uses. 

The idea of getting a strictly convex drawing by perturbing a convex drawing was pio- 
neered by Chrobak, Goodrich and Tamassia [Sj. They claimed to construct strictly convex 
embeddings on an O(n^) x O(n^) grid, without giving full details, however. This was improved 
to 0(n'^/3) X 0(n'^/3-) in pj. In this paper we further improve the "fine perturbation" step 
of jl2j to obtain a bound of O(n^) x O(n^) for grid drawings. Theorem ^ gives better bounds 
when the faces have few sides, and we allow grids of different aspect ratios (keeping the same 
total area). 

In the course of the proof, we need explicit (not just asymptotic) lower bounds on the 
number of primitive vectors in certain triangles. A primitive vector is an integer vector which 
is not a multiple of another integer vector; hence, primitive vectors can be used to characterize 
the directions of polygon edges. The existence of many short primitive vectors is the key to 
constructing strictly convex polygons with many sides. These lower bounds are derived in 
Section [3 based on elementary techniques from the geometry of numbers. 



2 



2 Preliminaries: Schnyder Embeddings of Three-Connected 
Plane Graphs 



Felsner fH| (see also Pllll|) has extended the straight-hne drawing algorithm of Schnyder, which 
works for triangulated planar graphs, to arbitrary three-connected graphs. It constructs a 
drawing with very special properties, beyond just having convex faces. These properties will 
be crucial for the perturbation step. 

Felsner's algorithm works roughly as follows. The edges of the graph are covered by three 
directed trees which are rooted at three selected vertices a, b, c on the boundary, forming a 
Schnyder wood. The three trees define for each vertex v three paths from v to the respective 
root, which partition the graph into three regions. Counting the faces in each region gives 
three numbers x, y, z which can be used as barycentric coordinates for the point v with respect 
to the points a, 6, and c. Selecting abc as an equilateral triangle of side length / — 1 (the 
number of interior faces of the graph) yields vertices which lie on a hexagonal grid formed by 
equilateral triangles of side length 1, see Figure^. Since / < 2n this yields a drawing on a 
grid of size 2n x 2n. 
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(a) (b) 

Figure 1: (a) A Schnyder embedding on a hexagonal grid and (b) on the refined grid after 
the initial (rough) perturbation 



This straight-line embedding has the following important property (see [HI Lemma 4 and 
Figure 11], [1 Fact 5]): 

The Three Wedges Property. Every vertex except the corners a,b,c has exactly 
one incident edge in each of the three closed 60° wedges shown in Figure [2h- 

From this it follows immediately that there can be no angle larger than 180°, and hence all 
faces are convex. Moreover, it follows that the interior faces F have the Enclosing Triangle 
Property, see Figure |^ ((SI proof of Lemma 7], |^ Lemma 2]): 

The Enclosing Triangle Property. Consider the line x = const through the point 
of F with maximum x-coordinate, and similarly for the other three coordinate 
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directions. These three hnes form a triangle Tp which encloses F. Then all 
vertices of F lie on the boundary of Tp, but F contains none of the vertices of Tp. 

It follows that interior faces with < 4 sides are already strictly convex. Throughout, we 
will call Tp the enclosing triangle of the face F. 

The Schnyder wood and the coordinates of the points can be calculated in linear time. 
Recently, Bonichon, Felsner, and Mosbah have improved the grid size to (n — 2) x (n — 2). 
However, the resulting drawing does not have the Three Wedges Property. An alternative 
algorithm for producing an embedding with a property similarly to the Enclosing Triangle 
Property is sketched in Chrobak, Goodrich and Tamassia 6 . It proceeds incrementally in the 
spirit of the algorithm of de Fraysseix, Pach and Pollack 7 and takes linear time. From the 
details given in ■6.1 it is not clear whether the embedding has also the Three Wedges Property, 
which we need for our algorithm. The original algorithm of Chrobak and Kant achieves 
a weak form of the Three Wedges Property, where F is permitted to contain vertices of Tp. 
Maybe, this algorithm can be modified to obtain the Three Wedges Property, at the expense 
of a constant-factor blow-up in the grid size. 




(a) (b) 

Figure 2: (a) Each closed shaded wedge contains exactly one edge incident to v. There may 
be additional edges in the interior of the white sectors, (b) A typical situation at a vertex 
which is perturbed. 




Figure 3: The three possible new positions for a single vertex in the rough perturbation. 
(Only the three boundary vertices a, b, c are pushed in directions opposite to these.) 



3 Rough Perturbation 

Before making all faces strictly convex, we perform an initial perturbation on a refined grid 
which is smaller by only a constant factor. This preparatory step will ensure that the subse- 
quent "fine perturbation" can treat each face independently. 
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(a) (b) (c) 

Figure 4: (a) A typical face F constructed by the convex embedding algorithm, (b) The new 
positions of the vertices of F which are pushed out are indicated, (c) The result of the rough 
perturbation. The perturbation of the vertices with question marks depends on the other 
faces incident to these vertices. 



We overlay a triangular grid which is scaled by a factor of 1 /7, see Figures 01 and El A 
point may be moved to one of the three possible positions shown in Figure 01 by a distance 
of \/3/7. The precise rules are as follows: A vertex v on an interior face F is moved if and 
only if the following two conditions hold. 

(i) The interior angle of -F at f is larger than 150° (including the possibility of a straight 
angle of 180°); and 

(ii) V is incident to an edge of F which lies on the enclosing triangle Tp- 

See Figure for a typical case. Such a vertex is then pushed "out" , perpendicular to the 
edge of Tp. We call the angle between the two edges incident to F and v the critical angle 
of V. For a boundary vertex different from a, c, the exterior angle is the critical angle, but 
these vertices are not subject to the rough perturbation. The three corners a, 6, and c are 
treated specially: they are pushed straight into the triangle by the rough perturbation, as 
illustrated in Figure ^ Examples can be seen in Figure jljj-c and Figure 01 The result of 




Figure 5: Example of the rough perturbation, 
perturbing the example in Figure^ is shown in Figure^. 
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(a) Case I 



(b) Case Ila 



(c) Case lib 



Figure 6: The cases in the proof of Lemma The figures show possible locations for the 
neighbors y and z oi x. 

There can be no conflict in applying the rules by regarding a vertex v as part of different 
faces: the bound of 150° on the angle, together with the Three Wedges Property ensures that 
there is at most one critical angle for every vertex (Figure [JJa). 

The result has the following properties: 

Lemma 1. After the rough perturbation, all faces are still convex. 

Moreover, if each vertex is additionally perturbed within a disk of radius 1/30, the only 
concave angle that might arise at a vertex v is the critical angle of v. 

Proof. It is evident that no critical angle can become bigger than 180°. For non-critical angles, 
this is also easy to see (cf. Figure Ht). (In fact, the second statement is a strengthening of 
this claim.) 

We now prove this second statement of the lemma by considering different cases. The 
reader who is satisfied with the existence of some small enough perturbation bound e > 
may skip the rest of the proof. We continue to show that we can choose £ = 1/30. 

Consider a non-critical angle yxz at a vertex x in a face F. We assume without loss of 
generality that x lies on the lower left edge i of the enclosing triangle Tp. 

Case I. The point x is incident to a critical angle of another face F' , and thus x is pushed 
out of F'. 

Without loss of generality, we can assume that x lies on the lower right edge of Tp/ , and thus 
X is perturbed in the lower right direction, as in FigureEK- (The other case, when x lies on the 
upper edge of Tp' and is pushed vertically upward, is symmetric.) By the definition of critical 
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angles, the angle in F' must be bigger than 150°. This excludes from F all points vertically 
above x or to the left of x. The upper neighbor z of x, which is a grid point, is therefore 
restricted to a closed halfplane right of a vertical line r at distance 1/2 from x. The lower 
neighbor y must lie on or above the line I that bounds the enclosing triangle Tp. Thus, y and 
z are restricted to the shaded area in Figure EK- Even if all three points are perturbed by the 
rough perturbation, they are still separated by a vertical strip of width 6?=^ — 2-^ = ^. 
An additional perturbation of ^ < 2^ cannot make the angle at x larger than 180°. 

Case II. The point x not perturbed by the initial perturbation. 

Case Ila. The point x has a neighbor on i. 
We can assume w.l.o.g. that it is the lower neighbor y, see Figure IHb- The angle yxz must be 
at most 150° because otherwise x would be critical. It means that z cannot lie to the left of x, 
and thus y and z are restricted to the shaded area in Figure IHh- Even if they are perturbed, 
they remain above the line s, which is obtained by offsetting the edge of the shaded region 
that is closest to x. The distance from x to s is 1/7 • y^3/7 0.0935 > Jj. Thus, there is 
enough space to additionally perturb the points x, y and z without creating a concave angle. 
(Actually, the vertex x will not even be perturbed in the fine perturbation.) 

Case lib. The point x has no neighbors on see Figure Et- 
This means that y and z lie on or beyond the next grid line parallel to i. The rough 
perturbation can move them closer to £, but they remain beyond another parallel line d." 
whose distance from x is 5/7 • y^3/4 ~ 0.618. This leaves plenty of space for additional 
perturbations of x, y, and z. □ 

After the rough perturbation, we will subject every vertex v that is incident to a critical 
angle to an additional small perturbation of a distance at most 1/30. The lemma ensures 
that, in order to achieve convexity at v without destroying convexity at another place, we 
only have to take care of one incident face when we decide the final perturbation of v. We 
can thus work on each face independently to make it strictly convex. 

4 Fine Perturbation 

We will now discuss how we go about achieving strict convexity of all faces. The rough 
perturbation helps us to reduce this task to the case of regularly spaced points on a line 
f Section 14. 1() . In Section we will describe in detail how the perturbed strictly convex 
chain is constructed for this special case. 

4.1 The Setting after the Rough Perturbation 

After the rough perturbation, we are in the following situation. Consider a maximal chain 
V2,V3, . . . , vk-1 of successive critical angles on a face F. These angles must be made strictly 
convex by perturbing them inside their little disks. (The two extreme angles at V2 and vk-i 
might already be convex.) The vertices ^2,^3,... ,vk-i he originally on a common edge of 
the enclosing triangle Tp, We first discuss the case when the vertices lie on the upper edge £ 
of Tp, forming a horizontal chain, as in Figure EK- (The extension to the other two cases is 
discussed in Section f4. 31 1 According to Lemma ^ we have to ensure that these critical angles 
are smaller than 180° after the perturbation. In Figure E^, these are the vertices V2, v^, and 
f4. Let us call these vertices critical vertices. In addition, we look at the two adjacent vertices 
vi and vk on F. By the choice of a maximal chain, they are not critical for F. They may lie 
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Figure 7: The setting of the fine perturbation process: (a) The initial situation after the rough 
perturbation. The angles in which it is necessary to ensure a convex angle are marked, (b) 
The circles in which the fine perturbation is performed. The size of the circles is exaggerated 
to make the perturbation more conspicuous, (c) A strictly convex polygon inside the circles, 
(d) The final result. 

on the same line as the critical vertices, as the vertices vi and in Figure EKj or they might 
lie below this line. To guide the perturbation of the points V2,. ■ ■ ,vk-i, we pretend that vi 
and vk are part of the chain, and we create surrogate positions v[ and v'j^ for these neighbors: 
First we move them from their original positions vertically upward to £; if they don't land on 
a grid point, we move them outward by 1/2 unit. Since the angles at V2 and vk~i are bigger 
than 150°, we are sure that v'i,V2, ■ ■ ■ ,vk-i,v'j^ lie on £ in this order. Finally, we subject v'l 
and v'j^ to the same rough perturbation as the critical vertices between them, and move them 
vertically upward. 

We place a disk of radius 1/30 around every perturbed point on this edge, including the 
two surrogate positions, see FigureEb. In the next step, to be described in Section we find 
a strictly convex chain which selects one vertex out of each little disk, as shown in Figure [7J;. 

This will make all angles at V2, ■ ■ ■ ,vk-i strictly convex. Finally, we use these perturbed 
positions for our critical vertices, but for vi and vk, we ignore their perturbed surrogate 
positions, see Figure I7|l. The true position of vi or vk may be determined by a different face 
in which it forms a critical angle (as is the case for in the example), or it might just keep 
its original position (like vi in the example). We only have to check that the angle at the 
left-most and right-most critical vertex {v2 and V4 in this case) remains convex: 

Lemma 2. Replacing the perturbed surrogate position v'l and v'^^ of the points vi and vk by 
their true positions does not destroy convexity at their neighbors V2 and vk-i in F . 

Proof. We first show that the rough perturbation does not actually perturb vi and vk to 
their surrogate positions v'^ or v'j^. It is conceivable that, say, vi lies on i and is perturbed 
upwards because of its critical angle in a different face F' , see Figure |H1 However, this would 
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contradict the Three Wedges Property for vi and F, creating two incident edges in a sector 
in which only a unique incident edge can exist. 
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Figure 8: A neighbor of a critical vertex cannot be perturbed in the same direction. 

Thus we conclude that vi and vk he below or on i, and they are either perturbed not at 
all or in a direction below £. 

Vertices V2 and V4 in the example of Figure [7| represent the possible extreme cases that have 
to be considered, f 5 represents a vertex that is pushed downward in the rough perturbation, 
and then subjected to a fine perturbation anywhere in its little circle. For visual clarity, the 
circles in FigureElhave been drawn with a much larger radius than 1/30. Since the circles are 
actually small enough, the angle at will be convex no matter where the point ^5 is placed 
in its own circle. (This position is determined when the critical face of ^5 is considered.) A 
similar statement holds at V2, where the perturbed surrogate position of vi in Figure [Tj: is 
replaced by the original position of vi; this will always turn the edge V2V1 counterclockwise 
and thus preserve convexity at V2. 

The argument works also for a chain of vertices on an exterior edge of the enclosing 
triangle. In this case, V2,V3, . . . ,vk-i are perturbed around their original position on i, 
whereas the neighbors vi and vk are moved inside the triangle and below i. Geometrically, 
the situation looks similar as for vertex vi in Figure except that vi is not pushed down 
straight but at a —30° angle. This movement is large enough to ensure convexity at ^2- □ 

4.2 Convex Chains in the Grid 

We have a number K of vertices = ai < 02 < • • • < ax < 2n — 1 on a horizontal line which 
form part of an array of 2n consecutive grid points. We want to perturb them into convex 
position. If the faces of the embedding have at most k sides, then K < k. It is more convenient 
to work with a rectangular grid. So we extend the hexagonal grid to a rectangular grid as 
shown in Figure IHl This grid will be refined sufficiently in order to allow a strictly convex 
chain to be drawn inside a sequence of circles. Figure gives a schematic picture of the 
situation. (This drawing is not to scale.) It is more convenient to discuss the construction of 
an upward convex chain. Inside each disk (of radius 1/30) we fit a square of side length 1/50, 
which is subdivided into a subgrid of width w and height h. More precisely, we are looking 
for a sequence of points pi = {xi,yi) in these circles, whose coordinates measure the distance 
from the lower left corner of the first circle in units of little grid cells. Two successive circle 
centers at distance 1 in terms of the original grid have a distance of S := 50w when measured 
in subgrid units. Thus we are looking for integer coordinates that satisfy a-i- S < x-i < ai-S+w 
and < Hi < h. Eventually, when the whole subgrid is scaled to the standard grid Z x Z, Xj 
and Hi will become true distances again. The total size of the resulting integer grid will be 
0{nw) X 0{nh). 
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(a) (b) (c) 

Figure 9: The hexagonal grid (a) is contained in a rectangular grid (b). A hexagonal grid twice 
refined (c) contains rectangular grids in three different directions. One of these rectangular 
grids is highlighted by thicker points. 

y 



i S = 50w 




w 



Figure 10: A convex chain formed by grid points in the circles. (Again, the radius of the 
circles is drawn much too large compared to their distance.) 



The convex chain pi,p2, ■ ■ ■ ,Pk has a descending part up to a point with minimum y- 
coordinate and an ascending part. We choose the two points with minimum y-coordinate to 
lie in the middle: We define M := [K/2\ + 1 and set yu-i = Um = 0. We will only describe 
the construction of the ascending chain from pM to the right. The left half is constructed 
symmetrically. 

The direction between two grid points is uniquely specified by a primitive vector, a vec- 
tor whose components are relatively prime. We now take a sequence of primitive vectors 
qi,q2, ■ ■ ■ , Qk-Mi Qi = {ui,Vi) with < Ui < w and Vi > 0, in order of increasing slope Vi/ui. 
Then we choose the difference vectors Ap as appropriate multiples of these vectors, in the 
following way. We have already defined yM '■= 0, and we choose xm arbitrarily within the 
permitted range of x-coordinates. Having defined pu+i-i, we define 

PM+i ■= PM+i-1 + S - Qi 

by adding as many copies of qi as are necessary to bring xm+i into the desired box: 

dM+i ■ S < XM+i < dM+i ■ S + W 

Since this box has width w, and Uj < w, this is always possible. 

We need K — M < K/2 primitive vectors qi (including the vector (1,0) from pu-i to 
Pm) The following theorem ensures that we can find these vectors in a triangle of sufficiently 
large area. 
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Theorem 2. The right triangle T = (0,0), {w,0), {w,t), where w > 1, w integer, and t>2, 
contains at least wt/A primitive vectors. 

The general proof is given in Section [51 We can however easily give an explicit solution 
for the special case t = 2 (corresponding to the choice w = k below, which leads to the most 
balanced grid dimensions): In this case, we can simply take the 1 + [u'/2j vectors (u),!), 
(«;-l,l), (k/21,1). 

We use Theorem |21 as follows. We choose an arbitrary width w < k for the boxes. By 
Theorem 121 we can set t := max{2,2K/w} to ensure that we find at least K/2 primitive 
vectors in the triangle T. The slope of these vectors is bounded by t/w. Let us estimate the 
necessary height h of the boxes. The last point px is connected to pM by a chain of vectors 
with slope at most t/w. The distance of x-coordinates is at most the width of the whole 
grid on which the graph is embedded, i. e., at most S ■ 2n = 0{wn); hence the difference in 
y-coordinates is at most t/w ■ 0{wn) = 0{tn) = 0{kn/w). It follows that the height h of the 
boxes is 0{kn/w). The total height of the resulting grid is 0{hn) = 0{kn'^ /w). 

This leads to part (i) of Theorem ^ Part (ii) is an easy corollary. As an extreme case, 
we can set w = 1 and perform only vertical perturbations. We get h < 2kn (without any 
additional constants depending on S) . 




(a) (b) (c) 

Figure 11: A rectangular grid (a), its 2 x 6 refinement (b), and a shearing (c) of the refined 
grid. Its grid-points coincide with the untransformed grid. 



4.3 Perturbation of Vertices on Diagonal Lines 

So far, we have treated only a sequence of vertices on a horizontal straight line. The same 
scheme can be applied to lines of the two other directions by applying the shearing transfor- 
mation (p 1-^ {y+^/2-x) (p ^ (j/-v^/2-z) '^hich moves points only in vertical direction. 
If /i is a multiple of w, the transformation will produce a grid like in Figure lllb which is 
contained in the original grid of Figure ITTb . For the range of parameters which is interesting 
for the theorem {w < k), the height h of the subgrid is never smaller than the width w; thus, 
the choice of /i as a multiple of w does not change the asymptotic analysis. One needs to 
reduce the size of the little square subgrid to ensure that the sheared square still fits inside 
the circle, and one has to adjust the quantity S accordingly. In addition, we have to select h 
and w as multiples of 14, to accommodate the grid of the rough perturbation and the refined 
rectangular grid of Figure O). All of this changes the analysis only by a constant factor. 

For the case of a uniform stretching of both dimensions {w = h), one referee has pointed 
out a simpler alternative method. After a blow-up by a factor of two, the original triangular 
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grid contains rectangular grids in all three grid directions, Figure IHt. Two further refine- 
ments by the factor 7 (for the rough perturbation) and then by the factor w are sufficient to 
accommodate the fine perturbation. 

On the exterior edges, the points must of course be perturbed to form an outward convex 
chain. 

For part (iii) of the theorem we have already mentioned that interior faces with k < 4 
sides are already strictly convex. If the outer face has 4 edges, it contains a single vertex on 
one of the sides of the outer triangle. The rough perturbation is thus sufficient to make the 
outer face strictly convex. 

The whole procedure, as described above, is quite explicit and can be carried out with a 
linear number of arithmetic operations. We calculate the 0{k) primitive vectors qi only once 
and store them in an array. Then, for every actual sequence of vertices on an edge, we can 
construct the perturbation very easily. The primitive vectors in the triangle (0,0), (^^,0), 
{w,t) according to Theorem [2 can be selected from the 0{wt) = 0{k) grid points in linear 
time with a sieve method. 

4.4 Numerical Experiments 

We have presented a general systematic solution for finding a convex chain by selecting grid- 
points from a sequence of boxes. One can find the optimal (i.e., longest) convex chain in 
polynomial time by dynamic programming, as described in more detail below. Results of 
some experiments are shown in the first column of Table ^ We restrict ourselves to the 
standard situation of selecting an n-gon from n adjacent boxes {K = n) which are squares 
{w = h). For several different sizes w, we computed the largest n such that a strictly convex 
n-gon can be found in a sequence of cells of size w x w. The factor (w + l)/n determines the 
necessary grid size w in terms of n. (By the convention of Figure [TUl a "tt; x w" grid consists of 
{w + 1)^ vertices; thus we give the fraction (w -|- l)/n instead of w/n.) Since the convex chain 
consists of a monotone decreasing and a monotone increasing part, connected by a horizontal 
segment in the middle, the necessary height w -|- 1 is at least 0.5 n. We see that this trivial 
lower bound is achieved for small values of n. The factor {w + l)/n increases with n, but not 
very fast. (The rectangular w x h boxes constructed in the proof of Theorem ^ would have 
w/n = 1, but h/n = 100.) 

The dynamic programming algorithm computes, for each point p in the w x h box, and 
for each possible previous point p' in the adjacent box to the left, the longest ascending and 
strictly convex chain (of length i) for which Pi-i = p' and pi = p. Knowing p' and p, it can 
be determined which points in the next box are candidate endpoints Pi+i of a chain of length 
i + 1. One can argue that, among these points Pi+i that are reachable as a continuation of 
p'p, only the w + 1 lowest points on each vertical line are candidates for endpoints Pi+i that 
form part of an optimal chain. Theoretically, the complexity of this algorithm is therefore 
0{w^h?). It turns out that, with few exceptions, every point p has only one predecessor point 
p' that must be considered: all other predecessor points pi-i have either a larger slope of 
the vector p — Pi^i or they are reached by a shorter chain. Therefore, the algorithm runs in 
0{wP'h) = 0{wkn) time, in practice. 

A simple greedy approach for selecting the points pi one by one gives already a very good 
solution: we choose Pi+i from the possible grid points in the appropriate box in such a way 
that the segment pj+i — pi has the slope as small as possible while still forming a convex 
angle at pi. The results in the right column of Table ^ indicate that this algorithm is quite 
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9.5227 
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11.2965 
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299,494 


13.3559 



Table 1: The length of the longest strictly convex n-gon in a sequence of square cells of size 
w X w, regularly spaced at distance S = bOw. 

competitive with the optimum solution. The running time is 0{kw). 

5 Grid Points in a Triangle 

In this section we prove Theorem |21 We denote by P := { (x,y) | gcd(x,y) = 1} the set of 
primitive vectors in the plane. 

It is known that the proportion of primitive vectors among the integer vectors in some 
large enough area is approximately 1/C(2) = 6/7r^ ^| Chapters 16-18]. Thus, a "large" 
triangle T should contain roughly 3/7r^ • wt ~ O.SOAwt primitive points. However, for very 
wide or very high triangles, the fraction of primitive vectors may be different. In fact, for 
t = 2, the bound wt/'i is tight except for an additive slack of at most 2. 

We will use special methods for counting primitive vectors when T is "very high" (i. e., w 
is fixed and below some threshold and t is unbounded. Section [5.1|) . when T is "very wide" 
(t is fixed and w is unbounded. Section l5.2() . and for the case when both t and w are large 
(Section l5.3j) . We use the help of the computer for the first two cases, but we use a general 
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bound for the last case. 



5.1 Fixed width, unbounded height 

For a fixed value of w, the function f(t) := |mP| can be analyzed explicitly. It is periodically 
ascending: 

f{t + w) = f{t) + C, 

where C = Y17=i '^(0 ^^^^ number of primitive vectors in the triangle (0,0), (w,0), {w,w), 
excluding the point (1,1). Euler's totient function denotes the number of integers 1 < 
j < i that are relatively prime to i, or equivalently, the number of primitive vectors (i,j) on 
the vertical line segment from (^,0) to — 1). 

The reason for the periodic behavior is that the unimodular shearing transformation 
(x, y) (x, y+x) maps the triangle (0, 0), {w, 0), {w, t), to the triangle (0, 0), {w, w), {w, w+t), 
which is equal to (0,0), {w,0), {w,w + t) minus the triangle (0,0), {w,0), (wjw). 

Therefore, it is sufficient to check that the "average slope" C/w of / is bigger than w/i, 
and to check 

fit) > tw/A (1) 

for the initial interval 2 < t < 2 + w. This can be done by computer: We sort all primitive 
vectors {x,y) with < x < w and < y/x < (w + 2)/'w hy their slope y/x. We gradually 
increase t from 2 to w + 2. The critical values of t for which must be checked explicitly 
are when a new primitive vector is just about to enter the triangle. 

We ran a lengthy computer check to establish for w = 1,2, . . . , 250 and for 2 <t < w+2 
(and hence for all t). In addition, we checked it for the range w = 251, 252, . . . , 800 and for 
2<t< 250. 



5.2 Large width 

In this section we prove Theorem |2l for small t and large w. T intersects each horizontal line 
y = i in a segment of length w — {w/t)i. In any set of i consecutive grid points on this line, 
there are precisely (/>(i) primitive vectors. We can subdivide the grid points on y = i into 
[{w — {w/t)i)/i\ > w/i — w/t — 1 groups of i consecutive points, leading to a total of at least 
{w/i — w/t — primitive vectors: 

l^nP|>i + E(f -y-i) ^« 

1=1 

For a given value of [t\ , one can evaluate the expression 

iTn PI > 1 + - f - 1),(,) > 1 + - ^ - (2) 

1=1 1=1 L ^ 

explicitly. The right-hand side of this bound is a linear function g{w): 
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Figure 12: (a) The triangle T' in Lemma IHl and its covering by squares, (b) The triangles T 
and T* (shaded) in Lemma 0J 

For example, for \t\ = 130, we have g{w) = w ■ 39.514 ... — 5153. It follows that g{w) > 
w ■ 131/4 > wt/A for w > 762. Performing this calculation by computer for [t\ = 6,7, ... , 130 
establishes Theorem [21 for 6 < t < 130 and w > 800. The interval 4 < f < 6 can be split into 
the ranges 4 < t < 4.5, 4.5 <t<5, 5<t< 5.5, and 5.5 < t < 6. For each range, we can use 
the above method with a tighter bound in ^ than t > [t\, and the estimate goes through in 
the same way. 

So let us consider the remaining interval 2 < t < 4: For 2 < t < 3, we can evaluate |rnP| 
explicitly: 

IT n P| = 1 + (u; + 1 - [f 1 ) + ([f 1 - [f - il ), (3) 

counting the primitive vectors on the lines y = 0, y = 1, and y = 2, respectively. For t > 3, 
the right-hand side of Q is still valid as a lower bound. We get 

|T n P| > 1 + + 1 - (f + 1)) + (f - (f - i + 1)) > - f ) 

The last expression is > wt/4: for 2 < t < 4. 

Thus we have proved the theorem for 2 < t < 130 and w > 800. 

5.3 Large triangles 

Lemma 3. Let T' be an axis-aligned right triangle of width a' and height b' , whose right angle 
lies on a grid point. Then 

areaT' <\T' r\I?\< areaT' + [a' J + [b' \ + 1 

Proof. This is simple. Suppose the right angle is at the right bottom corner of T' , see 
Figure 112b . Each lattice point in T' is the right bottom vertex of a unit square and these 
squares cover T' . To bound the area from below, we must subtract the squares which are not 
contained in T' . These squares form a monotone chain along the longest side of T' , and their 
number is \a'\ + [6'J +1. □ 



15 




Figure 13: If T* (shaded) contains no grid points, the triangle T' does not exist. 



Lemma 4. Let T he the right triangle (0,0), (a, 0), {a,b), with a,b > 1. Define T* as 
Tn{(x,y) : y > 1}. Then 

O'b , a ab , a 

a-b+—<\T*r\I?\<— + b-— (4 

2 25 - ' ' - 2 26 ^ ' 



In particular, 

nh 

<a + b 



\T*rM?\-- 



Proof. See Figure IT^ . The triangle T* has length a — a/b, height 6—1 and area ^(6 — l)(a — 
a/b) = ab/2 — a + a/ (26). Let T' denote the part of T* that lies left of the line x = [aj . This 
triangle contains the same grid points as T* . We assume first that T' is a nonempty triangle. 
The difference in areas lies in a rectangle strip of width < 1 and height 6—1: 

areaT* — (6 — 1) < areaT' < areaT* 

We can apply Lemma 01 to T' and obtain 

|r*nz2| = |T'nz2| < _ « + ^) + (a - ^) + (6 - 1) + 1, 

IT* n Z2| = |T' n I?\ > areaT' > (^^ - a + ^) - (6 - 1), 

from which the lemma follows. 

The triangle T' may not exist, as in Figure IT^ In this case, T*nZ^ = 0. Instead of arguing 
why the above derivation is valid also for this case, we establish the inequalities directly. Let 
6' > b — b/a denote the vertical extent of T at a; = [aJ. Then the fact that T' is empty is 
equivalent to 6' < 1. 

Then, from 1 > 6' > 6 — 6/a we conclude that a6 < a + 6. It follows that the lower bound 
in @ is at most 0: 

ab a , a + 6 a 

\- — - a - b < \ a-6<0 

2 26 -2 2 

The claimed upper bound in Q is always nonnegative, by the assumption 6 > 1. □ 

The number of primitive vectors can be estimated by an inclusion-exclusion formula, 
taking into account vectors which are multiples of single primes 2,3,5,7, .. ., vectors which 
are jointly multiples of two primes, of three primes, and so on, see |10l Chapters 16-18]: 

s s 
|TnP| = 1 + |T* nP| = i + ^fi{i) ■ |(i •T*)nz2| = i + ^;,(i) . \{lT)*nz^\ (5) 

i=l i=l 
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Here, fj,{i) is the Mobius function: = (—1)'^' if i is the product of k distinct primes and 
fj,{i) = otherwise. It is known that J^^^i = 1/C(2) = S/vr^, leading to the fact mentioned 
above that a fraction of approximately G/vr^ of the grid points in a large area are primitive 
vectors. 

Our sum in (jS)) goes to i = oo, but fov i > w or i > t, the set {jT)* n 7? is empty. 
Therefore, the formula is valid for S := mm{w, [t\}. We apply Lemma^and obtain 

1^ 1/1 wt ■ST^ Li(i) ^r^w + t 

|rnP| = i + j;M^).|(4rrnz2| >-^^- 

i=l i=l i=l 

wt f Q 1 

T 1^ " 5 



>-\-2-^]-Hs{w + t\ 



where Hg = 1 + 1/2 + 1/3 + -- - + 1/S is the harmonic number. The last inequality comes from 
bounding the remainder ^(^)/^^ — Yl'^s+i'^l^'^ ^ ^/^ °^ infinite series, whose 

value is 6/7r^. 

We distinguish the two cases for S: Case 1: w < t, and S = w. Then 

/ 6 1 \ . . / 3 1 2H„ 



2 \7r'^ w J yvr^ 2w w 

Case 2: w > t, and S = [t\. 

^, wt f Q 1 \ , 

|rnP|>-(^^-^j-ff,j(u. + <) 

Combining the two cases and setting n := min{t(;,t} gives 

3 1 2i/|^„j 



TnP > wt 



vr^ 2(n — 1) n 



Using the estimate Hi < j + ln(i + 1) with Euler's constant 7 ~ 0.57721, it can be checked 
that this factor is bigger than 1/4 for n > 250, thus proving the theorem for w,t> 250. 

On the other hand, the factor in ® is bigger than 1/4 for w > 800 and 130 < t < 250, 
proving the theorem also for this range. 



Wrap-up. The proof of Theorem [2 is now complete. On a high level, we distinguish three 
ranges for w: l<w< 250, 251 < -u; < 800, and w > 800. 

• Range 1: For 1 <w < 250, the theorem has been established in Section [5.11 

• Range 2: 251 < w < 800. For 251 < w < 800 and 1 < t < 250, the theorem has been 
established in Section l5?T] as well. For 251 <w< 800 and t > 250, it has been proved 
in Section [231 

• Range 3: Finally, for w > 800, there is a division into three cases: Section 15.21 takes 
care of the range 2 < t < 130. Section [^31 proves the bound separately for the ranges 
130 < t < 250 and t > 250. □ 
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6 Conclusion 



In practice, the algorithm behaves much better than indicated by the rough worst-case bounds 
that we have proved. We have not attempted to optimize the constants in the proof. For 
example, if we don't take a 7 x 7 subgrid but an 11 x 11 subgrid, and with a more specialized 
treatment of the outer face, the permissible amount of perturbation in Lemma H increases 
from 1/30 to 1/9, but it would make the pictures of the rough perturbation harder to draw. 

Bonichon, Felsner, and Mosbah ^ have used a technique of eliminating edges from the 
drawing that can later be inserted in order to reduce the necessary grid size for (non-strictly) 
convex drawings. This technique can also be applied in our case: remove interior edges as 
long as the graph remains three-connected. These edges can be easily reinserted in the end, 
after all faces are strictly convex. (For non-strictly convex drawings in the selection of 
removable edges and their reinsertion is actually a more complicated issue.) This technique 
might be useful in practice for reducing the grid size. 

Lower Bounds. The only known lower bound comes from the fact that a single convex 
n-gon on the integer grid needs r2(n^) area, see Barany and Tokushige or Acketa and 
Zunic ^121 for the easier case of a square grid. To achieve this area for an n-gon, one has to 
draw it in a quite round shape. In contrast, the faces that are produced in our algorithm have a 
very restricted shape: when viewed from a distance, the look like the triangles, quadrilaterals, 
pentagons, or hexagons of the n x n grid drawing from which they were derived. To reduce 
the area requirement below O(n^) one has to come up with a new approach that also produces 
faces with a "rounder" shape. 

Our bounds are however, optimal within the restricted class of algorithms that start with 
a Schnyder drawing or an arbitrary non-strictly convex drawing on an 0{n) x 0{n) grid and 
try to make it strictly convex by local perturbations only. Consider the case where n — 1 
vertices lie on the outer face, connected to a central vertex in the middle. The Schnyder 
drawing will place these vertices on the enclosing triangle, and at least n/3 vertices will lie 
on a common line. They have to be perturbed into convex position, as in Figures 171 or IIUI 

Let us focus on the standard situation when we want to perturb n equidistant vertices 
on a line, at distance 1 from each other. The n — 1 edge vectors pj+i — pi lie in a 2w x 2h 
box; they must be non-parallel, and in particular, they must be distinct. If Ay is the average 
absolute vertical increment of these vectors, it follows that Ay = Q,{n/w), and the total 
necessary height h of the boxes is r2(n(Ay)) = J7(n^/u;). Therefore, the total necessary area 
is Vl{hwn'^) = n{n'^). 

The argument can be extended to the case when only i7(n) selected grid vertices on a line 
of length 0{n) have to be perturbed. It can also be shown that our bounds in terms of k 
are optimal in this setting. The worst case occurs when there is a line of length n with n{k) 
consecutive grid points in the middle and two vertices at the extremes. 

Extensions. The class of three-connected graphs is not the most general class of graphs 
which allow strictly convex embeddings. The simplest example of this is a single cycle. 
A planar graph, with a specified face cycle C as the outer boundary, has a strictly convex 
embedding if and only if it is three-connected to the boundary, i. e., if every interior vertex (not 
on C) has three vertex-disjoint paths to the boundary cycle. Equivalently, the graph becomes 
three-connected after adding a new vertex and connecting it to every vertex of C. These 
graphs cannot be treated directly by our approach, since the Schnyder embedding method 
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of Felsner |H] does not apply. Partitioning the graph into three-connected components and 
putting them together at the end might work. 
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